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ABSTRACT 

The link between turbulence in star formatting environments and protostel- 
lar jets remains controversial. To explore issues of turbulence and fossil cav- 
ities driven by young stellar outflows we present a series of numerical simula- 
tions tracking the evolution of transient protostellar jets driven into a turbulent 
medium. Our simulations show both the effect of turbulence on outflow struc- 
tures and, conversely, the effect of outflows on the ambient turbulence. We 
demonstrate how turbulence will lead to strong modifications in jet morphology. 
More importantly, we demonstrate that individual transient outflows have the ca- 
pacity to re-energize decaying turbulence. Our simulations support a scenario in 
which the directed energy/momentum associated with cavities is randomized as 
the cavities are disrupted by dynamical instabilities seeded by the ambient turbu- 
lence. Consideration of the energy power spectra of the simulations reveals that 
the disruption of the cavities powers an energy cascade consistent with Burgers'- 
type turbulence and produces a driving scale-length associated with the cavity 
propagation length. We conclude that fossil cavities interacting either with a 
turbulent medium or with other cavities have the capacity to sustain or create 
turbulent flows in star forming environments. In the last section we contrast our 
work and its conclusions with previous studies which claim that jets can not be 
the source of turbulence. 

Subject headings: ISM: jets and outflows, ISM: clouds, stars: formation, turbu- 
lence, hydrodynamics 
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Introduction 



Both turbulence and outflows appear to be ubiquitous phenomena associated with star 
formation. Turbulence in molecular clouds is inferred from Larson's Laws which are e mpirical 
relat ionships between line-widths and size observed in many star forming regions ((Larson 
198ll ). Turbulence has also been inferred fro m direct measurement s of power spectra in 
molecular clouds across a wide range of scales (IHeyer &: Bruntl 12004 ) . Though the issue of 
whether the observed spectra are a conse quence of gravita tional fragmentation or mechanical 
forcing remains a topic of current debate (IField et al. Outflows driven by newly formed 

stars also appear uniformly across star forming environments. Outflows appear to form very 
early in the formation of a star (the Class phase) and continue until a sizable fraction of 
the stellar mass has been assembled. In addition, outflows occur across the entire spectrum 
of stellar masses fro m brown dwarfs (jCabritl 120071 ; IWhelan et al.l 120051 ) all the way to high 
mass O and B stars (IShepherdll2003l ). 



Given the apparent ubiquity of both turbulence and outflows, one can ask about their 
respective importance in terms of impact on the star formation process. In recent years 
the role of turbulence in moderating star formation h as become an active area of research 
( jMac Low fc Klessenl 12004 ; iKrumholz fc McKed 120051 ) . Some workers have argued that tur- 
bulent pressure in the form of a space filling isotrop ic distribution of eddies provides sup port 



for molecular clouds against gravitational collapse (jGoldreich fc Kwanlll974j ; lFleck!ll980l ). In 
this way turbulence is thought to be the agent which establishes equilibrium and provides a 
means for molecular clouds to be long lived structures. One difficultly with this view however 
is the recognition that both hydrodynamic and magnetohydrodynamic turbulence decays on 
a timescale of order the cloud free-fall time t Q ff (jStone et al.lll998l ; iMac Lowill999l ). Thus if 
clouds are stable, long-live d structures som e means must be found to support them against 
self-gravitational collapse ([Blitz et eil.l 120071 ) . 



When discussing issues of jet feedback one can distinguish between three scale lengths, 
"micro-scale" feedback concerns the effect of outflows on their own launch scale (L < 10 AU). 
"meso-scale" feedback relates to the impact of jets on environment associated with infalling 
envelopes (L < 10 4 AU). C onsiderable work has focused on these domains : i.e. the way 



outflows can sculpt envelopes dDelamarter et al 



ence on the accretion onto the star (IKrumholz et al. 



2000 



Cunningham et al.ll2005l ) or their influ- 



20051 ). Finally "macro-scale" feedback, 



which is the concern of this paper, is associated with the effect of outflows on the scale of 
star forming clusters or the entire molecular cloud itself. The role of outflows on large scales 
became of interest after it was discovered that some outflows extended to parsec scales which 
is of ord er cloud size or at least the size of the "clump" from which star clusters form withi n 
a cloud JBallv fc Devin"elll994j : iBallv et alill996l : Istanke et al1ll999l : iMcGroartv et al.l 120041 ). 
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Young star clusters such as NGC 1333 are the domains of interest in this problem. The 
impact of many stars, forming roughly co-evaly, on their parent cloud remains a subject of 
intense debate and this is preciously the location where issues in outflows and turbulence 
overlap. 

Numerous observational studies of star forming regions have shown that, in general, 
the energy in active outflows is large enough to s upport a cluster or cloud against collapse 



if that energy can be coupled to cloud tur b ulence (IBally et al.l Il996l ; iBally fc Reipurthl 12001 



Knee & Sandell 


2000; 


Quillen et al. 


2005: 


Warin et al. 


1996) 



the idea that outflows can provide a source of internal forcing to drive turbulence. Theoret- 
ical treatments of feedback from mul tiple spherical o u tflows creating a self-regulating star 
forming system was first explored by iNorman fc Silkl (Il980l ). Analytical work by Matzner 



flMatzner fc McKeel I2OO0I; iMatzner 



2001 



20071 ) have explored the role of collimated outflow 



feedback on clouds. IMatzner! (120071 ). in particular, developed a theory for outflow driven 
turbulence in which line-widths were predicted as functions of a global outflow momentum 
injection rate. iKrumholz et al.l (120061 ) have also considered the nature of feedback via out- 
flows, concluding that these systems provide an important source of internal driving in dense 
star forming cores. 



The first simulation based study of the problem was presented by iMac Low! (120001 ) 
whose results indicated that outflows could drive turbulent motions. Decay rates for outflow 
powered turbulence in this study were ~ 4x faster than that of uniformly driven turbu- 
lence. These early simulations we re of relatively low reso lutio n and further high reso lution 
work is needed. Recent work by iLi fc Nakamural (120061 ) and iNakamura fc Lil (120071 ) have 
mapped out the complex interplay between star formation outflows and turbulence. Their 
simulations include a self-consistent formulation of driving outflows from newly formed stars 
and they concluded that protostellar outflows were a viable means of generating t urbulence 
in star forming clusters. However, studies of single jets by iBanerjee et al.l (120071 ) came to 
the opposite conclusion. Exploring the volume filling averages of post shock material they 
concluded that single jets would not leave enough supersonic material in their wakes to act 
as a relevant source of internal forcing. 

Many studies of protostellar turbulence have focused on active outflows. These may 
not, however, be the only means of coupling outflows with turbulence. Numerical studies 
Cunningham et al.l (|2006al ) have shown that direct collisions of active outflows are not effec- 
tive at transferring jet/outflow momentum to cloud material (by active we mean that the 
momentum flux of the jet or wide angle win d rema ins roughly constant through out the in- 
teraction). Along similar lines iQuillen et al.l (120051 ) explored the relation between turbulent 
cloud motions and observed outflow activity in NGC 1333. Their results showed that the 
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location of active outflows did not effect the degree of turbulence in the cloud. iQuillen et al. 
( 120051 ) were, however, able to identify many sub-parsec sized, slowly expanding cavities. 
These they identified as dense wind-swept shells of gas that expand into the cloud after the 
driving source of the outflow has expired. The net mechanical energy required to open all 
these cavities represents a significant fraction of that required to drive the cloud turbulence. 
The authors argued that these cavities were relics of previously active molecular outflows. 
Because cavities have a longer lifetime than those opened by active and younger outflows, 
they would be more numerous than higher velocity, younger outflows. These cavities could 
provide the coupling between outflows and turbulence in molecular clouds. Thus one means 
of re-energizing turbulence in clouds may be through "fossil outflow" cavities which are dis- 
rupted af ter they have slowed to spe eds comparable to the turbulent velocity of the ambient 
cloud. In ICunningham et al.l (l2006bl ) high resolution AMR simulations of transiently driven 
outflows were carried out. These simulations, which included a treatment of molecular chem- 
istry and cooling, tracked the development of fossil outflow cavities driven by transient out- 
flows out to near parsec scales. The simulations provided observational predictions in terms 
of position- veloc i ty dia grams that were consistent with the momentum scaling relations used 



in 



Quillen et al.l (120051 ). 



While the potential for protostars to act as a means of internal forcing remains attractive, 
many outstanding questions remain open. The energetics of turbulent resupply must be 
understood in terms of a coupling coefficie nt between outflow and cl oud. The nature of 
outflow interactions must also be specified as lCunningham et al.l (j2006al ) have already shown 
that active outflows are not efficient it setting large volumes of gas into random motions. 
Thus the emphasis shifts to fossil cavities an their interactions with each other and/or with 
pre-existing turbulent medium. We expect that turbulence can be generated by collisions 
between fossil outflows because over the lifetime of a cluster the cavities will tend to fill 



the cluster volume I Cunningham et al.l (J2006aJ). Individual fossil cavities need not experience 
collisions but will still couple to the ambient medium. With regard to individual fossil 
outflow cavities and turbulence the mechanisms through their directed momentum/energy 
is "isotropized" ha s yet to be art i culate d. Something of a chicken and an egg problem exists 
in this regard. In I Quillen et al.l (120051 ) expanding outflow cavities where observed down to 
speeds a few times that of the turbulence v out fi ow ~ 2v tur b ~ 1 km s~ . They argued that 
the cavities would be subsumed by turbulent eddies once they decelerated to speeds of order 
v turb- The expectation that the directed momentum and energy of a single decelerating cavity 
can be given up to turbulence must be explored. This is the question we seek to address in 
this paper. 



In the current work w e perform a series of a simulations meant to take the work of 
Cunningham et al.l (l2006bl ) to the next level of complexity. This paper is meant as an critical 
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intermediary step to full scale simulations of multiple interacting transient outflows. Here 
we consider the role of individual transient jets which are ejected into a decaying turbulent 
medium. We note that the evolution of a radiative jet in a turbulent medium has yet to be 
explored and our simulations are relevant to YSO jet morphology and the reenergization of 
turbulence via jets and the disruption of the outflow cav i ties. Finally our simulations are 
relevant to issues addressed in the work of iBanerjee et al.l (120071 ) and their conclusions that 
jets can not drive supersonic turbulence. We note that the question of the steady state rate 
of decay of such flows remains outside of the scope of the present work as this question could 
only be addressed by models that include a steady injection of multiple outflows which are 
evolved for longer timescales than the nu merical models consi dered in this paper. This issue 
will be addressed in a forthcoming paper ([Carroll et al.ll2008l ). 



2. Numerical Model 

We have carried forward a series of numerical simulations to investigate the interaction of 
young stellar jets with a turbulent environment. Each of the simulations is carried forward 
on a periodic domain spanning L x x L y x L z = 4000 x 2000 x 2000 AU discretized to 
a 504 x 252 x 252 computational grid. An environment of density p a = 50 cm -3 with 
temperature T a = 200 K is set into a state of decaying turbulence from an unspecefied source. 
Possible candidates for supplying background turbulence within star forming environments 
include an initial turbulent state present du ring the formatio n of the cloud from compressive 
motions in a turbulent interstellar medium (lElmegreenll2000l ) , differential rotation of galactic 
disks, gravitational collaps e, the ionizing ra diation, winds and supernovae from massive starts 
and protostellar outflows (IMac Lowlll999l ). The spectrum of initial motions can be written 



as: 



with wave numbers 



H,j,k 



y^-j ft A x sin(kj j^ • x + <p x ij^k) 
Sj,j,fc A y ij,k sin(k i):; - jfc ■ y + 4> y i,j,k) 
52i,j,k Az i,j,k sin(k i j i jfe ■ z + <j) z iij:k ) 



2ni 



x + 



2nj 



y + 



2nk 



z. 



2000 AU 2000 AU " 2000 AU 

The summation is taken over all i, j, and k between and 31, excluding the % = j = k = 
term. The longest wave-mode seeded in each spatial direction is therefore equal to the length 
of the short edges of the computational grid and the shortest wave-mode in each spatial 
direction is resolved across approximately 16 computational zones. The phase angles <f>*,i,j k 



- 6- 



are chosen at random between and 2ir. The direction of each wave-mode is determined 
by choosing two of the amplitude components, A x and A y at random between —1/2 
and 1/2. The third amplitude A z is computed to satisfy the constraint that the initial 
velocity field be solenoidal. Each velocity amplitude A^jk is normalized so that the initial 
velocity power spectrum follow a power-law dependence on wavenumber £?jfo(|k|) ~ |k| _/3 
with net kinetic energy prescribed by the desired root mean squared (RMS) turbulent speed 
(vturb) as V J E Ko (k)dk = | PaVturb where V is the volume of the computational domain, 
Ek()z) is the power spectrum of velocity perturbations, 

E K (k) = |^K(x))| 2 + |^K(x))| 2 + |^(i; 2 (x))| 2 , 

and T is the Fourier transform operator, = j / x e lkx c?x. It is also convenient to define 
the one dimensional power spectrum as 



E K (k) = J E K {k)5{\k\ - k)dk. 



The simulations which begin with an initially turbulent flow field are all seeded with an initial 
RMS turbulent velocity v tur b = 10 km s _1 , yielding an initial RMS turbulent Mach number 
of 6.03, and velocity pow er spectrum index of = —2, cha racteristic of highly supersonic 



Burger 's-type turbulence (IPorter et al.lll992l ; Mac Lowll2003l ). 



Our simulations employ periodic boundary conditions at each boundary interface. This 
choice of boundary conditions prevents boundary artifacts from polluting the spectrum of 
motions in the grid. Jet-driven cavities are therefore launched from a grid-embedded region 
rather than from grid boundaries. This choice of domain setup allows the simultaneous use 
of periodic grid edges and jet launching. The jet launch region extends from 125 AU < 
x < 250 AU, and (y - 1000 AU) 2 + (z - 1000 AU) 2 < r) where rj = 250 AU. The launch 
region maintains a constant density pj = 250 cm" 3 and temperature Tj = 200 K profile and 
enforces a time dependent jet velocity as 



v x (t) 



v D exp {-t/t ff) t < 2t ff 
otherwise 



with initial velocity v Q = 150 km s~ . The jet velocity is smoothed following a quadratic 
form so that the outer radius of the launch region have a velocity that is 90% of the velocity 
at the center of the jet. 

A total of five simulations have been performed using the AstroBEAR astrophysical fluid 
code to carry forward the integration of the Euler equations using t he MUSCL-Hancoc k shock 



capturing scheme and Marquina flux function as described by ICunningham et al.l (120071 ) 



with an energy sink term which models the effects of optically thin radiative energy loss 
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as described by ([Cunningham et al.l l2006al ) . These simulations did not use the AMR or 
magnetic field capacities of AstroBEAR. The suite of simulations include three models in 
which jets of with varying decay parameters are driven into a decaying turbulent environment 
(Run 1-2), one model which follows the decay of the turbulent medium in the absence of 
disruption by a jet (Run 0) and one model which in which a non-decaying jet is driven into 
a quiescent environment (Run J). In table HJ we summarize the initial conditions and results 
each simulation. The fifth column from the left tabulates the cooling efficiency as measured 
by the net energy loss at the end of the simulations, the sixth column tabulates the total 
energy injected into the grid at the jet launch region over the duration of the simulations 
relative to the initial turbulent energy seeded in the grid and the rightmost column tabulates 
the least squared fit to the velocity spectrum power law index at the end of the simulation 
P\t=tf for 10 15 AU < k < 10 2 3 AU. In our simulations the power-law spectral dependence 
expected for inviscid turbulent flow begins to break for larger k due to numerical dissipation 
and for smaller k due to the inability to represent modes with characteristic length greater 
than the shorter edges of the simulation domain. 



Table 1: Simulation Summary. 





f off 




Turbulence 


Final Time tf 


Cooling Efficiency 




P\t=t f 


Run 


no jet 




/ 


922 yr 


36.1% 




-2.09 


Run 1 


8 yr 




/ 


922 yr 


49.0% 


0.285 


-2.24 


Run 2 


100 yr 




/ 


221 yr 


9.45% 


3.56 


-2.38 


Run 3 


non- decaying 


jet 


/ 


161 yr 


22.7% 


71.0 


-2.43 


Run J 


non- decaying 


jet 




161 yr 


22.1% 




-2.46 



The simulations presented in this work are of considerably shorter timescales than that 
of actual protostellar outflows. This is primarily due to computational constraints imposed 
by the need to simultaneously resolve the overall outflow, the outflow injection radius Rj 
and turbulent eddies with scale length < Rj. The simulations of extinct outflow sources 
presented here are interpreted as models of the fossil relics of extinct protostellar outflow 
activity in molecular clouds. The simulated of fossil cavities are intended to be roughly 
comparable to actual protostellar outflow cavities at comparable time relative to the outflow 
shut-off time, t/t ff. 



3. Results 



The right panel of figure [T] shows a crosscut about the mid-plane of the simulation 
domain at the end (t = tf) of each simulation. The initial and jet launch conditions of each 
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simulation are summarized in the left three columns of table [TJ The right three columns of 
table [H list diagnostic parameters at the end of each simulation. The "cooling efficiency" 
which is taken as the fraction of the total energy budget lost from the simulation domain via 
radiative cooling is listed for each run. This quantity measures the radiative loss of energy 
from the system that may have otherwise been available to sustain turbulent motion. Note 
that this cooling efficiency increases for simulations that run for longer physical time. This is 
because these simulations have more time for radiative energy loss to occur. It also increases 
for simulations with longer jet decay times as these jets feed energy into their bow shocks 
for longer times and have stronger cooling compared with jets that turn off earlier. 

The propagation speed of the Mach disk structure at the head of the simulated outflows 
can be estimated anal y tically for the case of the continuously driven jet models (run 3 and 



run J). iBlondin et al.l (119901 ) consider the momentum flux across the inward and outward 
facing surfaces of the Mach disc structure at the head of jet-driven flows to derive an analytic 
estimate for jet bow-shock propagation speeds, 



1 + 



EL 

Pa 



This expression evaluates to Vb s = 104 km s -1 for the continuously driven jet model param- 
eters considered in this paper. The shocks delineating the Mach disk at the head of the 
continuously driven jets (run 3 and run J) propagate ~ 3500 AU over tf — 161 yr, cor- 
responding to propagation speeds of 103 km s _1 , in excellent agreement with the analytic 
estimate. Comparison of the two simulations shows that the turbulence does not effect the 
average propagation speed of driven jets. 

Figure [2] shows the time evolution of the total mechanical energy in the grid for each 
simulation. Because the mechanical energy of the continuously driven jet quickly overwhelms 
that of the turbulent environment, the time evolution of the mechanical energy for run 3 and 
run J are indistinguishable in figure [2j For decaying jets, run 1 and run 2, the net mechanical 
energy increases until the driving source begins to decay at t ~ t ff- Interestingly, for t > t ff 
the mechanical energy decays in a roughly similar form as that of the control simulation, 
run 0, which which follows the decay of the turbulent environment alone 
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Fig. 1. — Velocity power spectra at several epochs (left) and crosscuts about the mid-plane 
showing the logarithm of the gas density, log(p[cm~ 3 ]), at the end of the each simulation. 
The panels show RunO, Runl, Run2, Run3, and RunJ from top to bottom. 
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Fig. 2. — Evolution of the total kinetic energy in the computational domain for each sim- 
ulation (left) and the evolution of the kinetic energy for runO and runl only (right). The 
kinetic energy of both run J and run 3 are dominated by the continuously driven jet and the 
two curves are indistinguishable in the left plot. Note the similar form of decay for the jet 
with turbulence cases (run 1 and run 2) to that of the turbulence-only control (run 0). 



3.1. Jet Morphology 



The morphological effects of the turbulent medium on the jet-driven bow-shock in run 
3 is readily apparent despite the comparatively slow turbulent speed in the ambient medium 
(vturb is only 10% of the jet propagation speed). This is particularly apparent in the change 
in behavior of the vortex ring in run J located near x ~ 2750 AU in comparison to run 3. 
Specifically the formation of a clearly defined vortex ring is suppressed for the case of the 
identical jet conditions which are driven into a turbulent medium in run 3. Furthermore, the 
entire bow shock from the base of the jet to its head shows significant disruption compared 
to the control case. This is most apparent, and can be most easily analyzed quantitatively, 
at the Mach disk at the jet head. In the jet which is driven into a quiescent environment 
(run J) a thin shock-bounded Mach disk structure forms at the head of the outflow. In 
the case where the same jet launch conditions are imposed into a turbulent environment 
(run 3), the Mach disk is observed to broaden significantly. Figure [3] presents a comparison 
of the Mach disk width (Wm) f° r both runs. The figure was constructed by measuring the 
extent of the Mach disk at its widest point using compression in the jet head as a tracer 
(p > 750 cm -3 ). The figure shows the Mach disk width increasing in the turbulent simulation 
by more than a factor approximately 4 over that in the control case. 
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To quantify this change run J shows a Mach disk width of Wm ~ 250AU. Dividing 
by the age of the jet yields an average expansion rate of the Mach disk of v exp = Wu/tf = 
7.1 km s~ . For run 3, the jet driven into the turbulent environment, the Mach disk structure 
shows Wm ~ 800AU. The average expansion rate is v exp = W M /t f = 22.7 km s" 1 . Note that 
the additional expansion of the Mach disk in the turbulent case exceeds the initial ambient 
turbulent velocity, v tur b = 10 km s~ . Thus it appears that more than just the entrainment 
of ambient turbulent motion is at play in this behavior. 

The change in the behavior of Mach disk highlights the fact that even structures which 
are denser and propagate more quickly than the ambient turbulence are significantly dis- 
rupted in its presence. This result suggests an interpretation of jet evolution in a turbulent 
medium which can account for the eventual disruption of fossil cavities in decaying jets. It 
is well known that initially laminar shock flows can be di srupted via the action of various 
instabilities associated with strong c ooling: i.e . radi ative ([Sutherland et al.l 120031 ) and thin 
shell modes IVishniac fc Ryul (119891 ); IVishniad (11994J ) . Given this susceptibility to fluid in- 
stability, the motion of the ambient turbulence with eddies occurring at a variety of scales 
should provide a space filling environment of multi-mode perturbation seeds. Because the 
turbulent field is composed of a continuous spectrum of eddies it will seed the outflow shell 
with perturbations which correspond to the fastest growing modes of any instability which 
the outflow-swept shell may be subject. These seeds act to initiate instability growth in the 
bow shock delineating the cavity and facilitate the coupling of the energy and momentum 
of the original outflow to the turbulent environment. 

Consider for example a shell of shocked gas of width h with internal sound speed c. The 
growth rate for modes o f the Non-Linear T h in Sh ell instability for an initial displacement L 
can be approximated as iHueckstaedt et al.l (120061 ): 



CJ 1/2 ckVkL 



Where is of order 1. The growth of the instability requires the displacement L be of 
order the shell width L ~ h. Since the shells are thin h << rj and the time scale for the 
bow-shock to fragment will be short compared to a characteristic hydrodynamic timescale of 
the flow which in this case we could take to be the crossing time of the jet across the cluster 
scale t c i ~ r c i/vj. Thus considering only the non-linear thin shell instability with L ~ h and 
k ~ 1/h we see that 



c r c [ 

rt c / ~ -r ~ 

vj h 



M~ 



h ' 



(2) 



Given typical Mach numbers in YSO jets of M ~ 100 and r c i ~ 1 pc we have Tt c i » 1. 
Thus perturbation seeds provided by the turbulence will lead to the fragmentation of the 
bow shock by the time the cavity has reached is maximum extent. 




t [yr] 



Fig. 3. — Temporal evolution of the Mach Disk height of the continuously driven jet in 
a turbulent environment, run 3 (solid line) and the continuously driven jet in a quiescent 
environment, run J (dashed line) plotted on a linear-linear scale. 



Putting the issue of the generation of turbulence aside for a moment we see that issues 
associated with the propagation of active jets are raised when we consider the disruption 
of the bow shock. Figure H] shows a three dimensional visualization of run 3. Note the 
fragmentation of the outflow into clumpy filaments. In particular, note the fragmentation of 
the Mach disk associated with the jet bow-shock complex into clumps. Each clump drives 
a smaller scale bow-shock in its wake, resulting in the formation of filamentary density en- 
hancements in the Mach disk. Similar irregular structures, including clumps and filaments 
in the Mach disk, have been ob served in a numb er of jet systems inc l uding HH 1/2 and 
HH 47 as h a s bee n discussed by iBally et al.l (120021 ) and lHartigan et al.l (120051 ) respectively. 
Yirak et al.l (120071 ) have noted similar bow-shock structures when the head of a propagating 
jet is impeded by stationary density homogeneities embedded within the ambient environ- 
ment. Thus it appears that filamentary fragmentation of protostellar jet bow-shocks may be 
characteristic of an inhomogeneous environment in general, regardless of the form of such 
inhomogeneities. 



Observations by IBally &: Reipurthl (120031 ) and IBally et al.l (120061 ) also reveal shock de- 
lineated knotty deflected structures within protostellar jet beams. Such deflections are puz- 
zling when they appear disc ontinuous, with th e jet showing a sharp deflection only some 
ways down the beam. While I Yirak et al.l (120071 ) find that sufficiently dense inhomogeneities 
can cause similar jet deflections, the motions of the turbulent environment in the present 
simulations show no comparable effect on the continuously (run 3) or long-driven (run 2) 
jets. However, the turbulent motions of the ambient environment do sufficiently disrupt the 
lateral edges of the jet bow-shock structure to expose the driving jet beam to the influence 
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of nearby turbulent eddies. It is interesting to note that the turbulent eddy scales which 
are represented in the computational domain (L < 4000 AU) do act to compress the edges 
of the long decay time jet beam (run 2). Thus it is possible that supra-jet-scale flow eddies 
(L ^> bow-shock radius) which would have higher velocity may drive bending and global 
disruption of jet stru ctures. On the issue of large-scale disruption, We note the work of 



Lebedev et al.l (120041 ) in which the deflection of supersonic jets via interaction with a large- 



scale transverse flow was considered. This study explored the issue both experimentally and 
numerically and addressed the issue of a jet interacting with side-wind after propagating 
some distance from the jet source. These studies raise the possibility of an encounter with a 
large eddy as the cause of discontinuous jet bending seen in observations. However, the effect 
of very large scale eddies are not considered in the present study as their inclusion would 
require a much larger computational volume and many more grid points to simultaneously 
achieve adequate resolution of the jet and the large scale flow. 



3.2. Turbulence and Kinetic Energy Power Spectra 

The left column of figure [1] shows plots of the velocity power spectra for each simulation 
at several times over the course of the simulations. The power spectra have been constructed 
by binning the three dimensional Fourier transform of the computational domain into spher- 
ical slices of spectral width Ak = 2n/Ax = 0.792 AU -1 . 

First we note the degree to which the outflows interact with the background turbulence. 
Only the jet driven cavity with the shortest driving decay time in run 1 becomes completely 
disrupted by, and subsumed into, the background turbulence. However, the more slowly 
propagating radial edges of the cavity bow-shock in the longer jet decay time runs all show 
some degree of coupling with the ambient turbulence. In the case of the continuously driven 
jet in the turbulent environment (run 3), approximately one cavity bow-shock radius (Rbs) 
crossing time (t BS ~ Rsshturb) has elapsed by the end of the calculation. In that simulation 
the oldest portion of the outflow limb, which is closest to the outflow source, shows the most 
disruption qualitatively. 

Moving now to the power-spectra we see the longer-lived jet sources of run 2 and run 3 
show increasing power in the largest flow eddies (smallest wavenumber). This indicates that 
the outflow cavities in these simulations provide sufficient power at length scales comparable 
to the width of the simulation domain to support turbulent motions in the ambient flow 
against energy decay. Even though the rapidly decaying jet (run 1) injects comparatively 
less energy, we also find more energy in the largest eddies compared to the case of turbulence 
without a jet. Much of the mechanical energy injected into the domain by the outflow source 
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Fig. 4. — Three dimensional volume rendered realization of the logarithm of the gas density, 
for the continuously driven jet in a turbulent environment with semi-transparent isosurfaces 
rendered at p = 2.33, 3.67, and 316 cm -3 . 

appears as large scale motion (smallest wavenumber) in this simulation as well. 

The increasing power at small wavenumber exhibited by the jet models can be inter- 
preted as the net result of the decay of turbulent energy at large wavenumber along with 
the resupply of turbulent energy by the injected outflow toward smaller wavenumber. The 
more powerful outflow source in run 3 resupplies energy toward small wavenumbers faster 
than the turbulent cascade of drains energy toward larger wavenumbers. In the case of the 
outflow cavity models (run 1 and run 2), modes with the smallest wavenumber increase even 
as the total kinetic energy in the grid decreases (figure [2]). We therefore identify the driving 
scale of turbulent resupply associated with the scale of cavity propagation as k drive ~ 1/L BS . 
Furthermore, the simulations that drive outflows exhibit steeper power-law slopes than the 
decaying turbulence case. At the end of the respective simulations, the decaying turbulence 
model in run reveals a spectrum with power law index f3 = —2.09 compared to the jet-only 
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model with j3 = —2.46. These cases bracket the slope of the spectra for the outflow with am- 
bient turbulence models which show a clear steepening with increasing outflow energy (Table 
[TJ). This steepening can be interpreted as an outflow-induced suppression of sub-cavity scale 
flow eddies. The power spectra therefore indicate that outflow cavities act to support the 
ambient turbulence in two ways. First the outflows power eddies of comparable extent to 
the length of the cavity Lbs and second, the opening of outflow cavities inhibit the cascade 
of energy to sub-outflow scales. 

We now address the issue of cavity disruption. Note that density distribution of run 3 in 
the neighborhood of the jet is characteristic of a uniform, laminar flow. After a short time, 
the energy injected by the jet greatly exceeds that of the ambient turbulence. Therefore, 
the power spectra at the end of this simulation is dominated by the jet at all scales and the 
spectra of run 3 appears indistinguishable from that of the jet-only control in run J. The 
density distribution crosscuts shown in the right column of figure [lj however, reveal that 
the outflow cavity driven by the short, nearly impulsive jet in run 1 has been completely 
subsumed into turbulent eddies. By the end of the simulation the flow field for the short jet 
pulse (run 1) is qualitatively simil ar to that of the turb ulen ce-only control (run 0). Th is is 



the behavior that was expected by iQuillen et al.l (120051 ) and ICunningham et al.l (j2006b[ ) for 



long extinct outflow structures embedded in turbulence. Furthermore, the net mechanical 
energy of the disrupted outflow cavity in run 1 decays at a rate comparable to that of the 
turbulence-only control in run 0. This implies that the disrupted outflow cavity becomes 
turbulent itself as it evolves under the influence of the turbulent environment. 

The time dependent evolution of the outflow cavity driven by the short jet pulse is 
illustrated in figure [5j The figure is composed of several plots which show cross cuts about 
the simulation mid-plane of the logarithm of the Mach number (left), the logarithm of the gas 
density in cm -3 (center), and the spatial distribution of the fraction of outflow ejected gas as 
followed by an advected contaminant (right) at t — 92, 230, 461, 691, and 922 yr from top 
to bottom. The color mapping of the plots of Mach number is constructed to reveal regions of 
supersonic flow only. The images in the right column showing the advected "jet contaminant" 
reveal that by the end of the simulation the ejected gas has been widely distributed by the 
turbulent flow across the computational domain. The end of the simulation at t — 922 yr 
corresponds to approximately one turbulent crossing time across the shorter edges of the 
computational domain of L/v tur b ~ 200 AU/10 km s _1 = 948 yr. By approximately one 
bow-shock radius crossing time, t ~ 461 to 691 yr, the initially laminar jet flow has evolved 
into a turbulent flow indistinguishable from the turbulent motions of the cloud. Furthermore, 
the flows with the highest Mach number at the end of the simulation correspond to regions 
containing the highest concentration of jet gas. This indicates that one crossing time after 
the outflow source expires the transform of bulk, directed fossil cavity motion into a fully 
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turbulent flow is complete. 




Fig. 5. — Evolution of the short-lived jet simulation, run 1 showing cross cuts about the 
simulation domain mid-plane of the logarithm of the Mach number (left), the logarithm of 
the gas density in cm -3 , (center) and the spatial distribution of the fraction of outflow ejected 
gas as followed by an advected contaminant (right) at t — 92, 230, 461, 691, and 922 yr 
from top to bottom. By the end of the simulation the jet-launched gas is widely distributed 
throughout the domain, the outflow cavity has been thoroughly subsumed into the ambient 
turbulence and much of the domain retains supersonic flow speed. Furthermore, the flows 
with the highest Mach number at the end of the simulation correspond to regions containing 
the highest concentration of jet gas. 
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3.3. Protostellar Jets and the Production of Turbulence 



The question of if, where and how protostellar outflows drive turbulence in their en- 
vironments remains an open issue. I n this work we have a ttempted to address one facet 
of the problem. In the recent work of iBanerjee et all (120071 ) the problem of outflow driven 
turbulence was also addressed however these authors came t o a ve ry different conclusion 
then we have in this and previous works. In IBanerjee et all (120071 ) single outflow cavity 
structures where used as a means of determining the efficacy of outflows as a source of su- 
personic turbulence. In that work it was argued that regardless of the momentum imparted 
to the cloud by an outflow, effective outflow-turbulence coupling will only occur if a sig- 
nificant volume fraction of the outflow cavity retains su personic flow speed at the time it 
is disrupted into the cloud. Using numerical simulations IBanerjee et all (120071 ) found that 
most of the volume overrun by either transient or continuously driven isothermal jet models 
remain subsonic after passing trough the bow-shock and thereby conclude that outflows are 
unlikely to act as a significant source of supersonic turbulence in molecular clouds. The 
results in our paper point towards the opposite conclusion. In what follows we discuss the 



differences between the basic scenario for jet driven turbulence presented in Banerjee et al. 



(2007|) a nd the ideas which u n derpin this work an d our pr evious studies ICunningham et al. 



(j2006al ); iQuillen et all (120051 ); ICunningham et all (l2006bT) . In addition we address specific 
issues raised by the simulations in Banerjee et al ( 20071 ) 



The basic premise of our work is that fossil cavities driven by transient protostellar 
outflows couple strongly to their surrounding cloud environments and drive turbulence within 
those clouds. We conjecture that the coupling can occur in two ways. First to produce 
turbulence in an initially quiescent medium interactions in the form of collisions between 
fossil cavities are required. Second, once fully developed turbulence exists individual fossil 
outflows can still give their energy and momenta back to the turbulent flow through the 
disruption of the cavity shell. The simulations presented here address the second process. 
With regard to the initial gener ation of turbulence via inte racting fossil cavities we note a few 
key points. First, as shown in ICunningham et all (j2006al ). in a typical cluster environment 
one can expect that every parcel of ambient gas will be overrun by an outflow at least once. 
One can also estimate occurrence of collisions. For clouds containing a volumetric density 
of iV* protostellar driven bipolar outflows each of radius R, length L and lifetime t out the 
probability P of cavity overlap increases as: 



p = nSr 2 l- 



t 



out 



tel 



oud 



(3) 



where t c i ou d is the age of the cloud which we assume fixed. Th us the older the outflow the 
higher the probability that cavities will overlap. In the work of ICunningham et all (j2006al ) 
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direct numerical simulations showed that active outflows, where jet beams collided, were 
less efficient at stirring ambient material than the sum of two similar outflows which do not 
collide. This was because the beams could not pass through each other and forced a shock 
mediated redirection in which significant kinetic energy of the flow was lost to radiation. 
The conclusion of that study was not, however, that collisions of all kinds where ineffective. 
Instead the authors conjectured that interactions between slower, longer lived (and larger) 
fossil cavities could be the source of turbulent driving. Indeed the important difference 
between active outflows and fossil outflow cavities is in the later case overlapping outflows 
can penetrate through each other causing disruption but not redirection. This is due to the 
fact that the fossil cavities are hollow and contain no active momentum bearing jet beam. 
We note that the work of IMatznerl (120071) provide s an a nalytical description of just such 
a process. Using dimensional arguments IMatznerl (120071 ) derives a length scale at which 
impulsively driven (spherical) outflows will interact. Given an ambient cloud of density p 
and volumetric outflow rate S and with each outflow driven with momenta I Matzner finds 
turbulence will be driven at a scale length: 



J1/7 



J turb 



^1/751/7 



(4) 



Typical values of p, I and S appr opriate for NGC 1333 y ield a n interaction sca le L tur h ~ 1 pc. 
Simulation studies presented in iLi fc Nakamural (120061 ) and iNakamura fc Lil ((2007J) support 
the basic scenario tha t mult iple outflows can drive turbulence in a protstellar cloud. In 
addition ICarroll et al.l (120081 ) have directly confirmed Matzner's prediction with a a longer 
effective scale length obtained due to jet collimation of L e ff w 3L turb . Thus we conclude that 
multiple overlapping, interacting fossil cavities can drive turbulence in an initially quiescent 
medium. 

In this paper we have shown that interactions need not occur for a fossil cavity, already 
in a turbulent medium, to resupply its environment with turbulent energy. Our identification 
of the disruption and eventual subsumption of the fossil cavities demonstrates the pathways 
through which an outflow interacting with only ambient turbulence can act to support a 
turbulent cascade. This can been seen in the time evolution of the Mach number of the short 
jet pulse driven outflow cavity (run 1) shown in the left column of figure revealing that the 
outflow cavity becomes significantly disrupted and subsumed into the ambient turbulence 
before becoming subsonic. Thus in the present simulations, the seeding of flow instability 
via ambient turbulence provides the mechanism for disruption. 

The scenario through which prot ostellar outflows drive turbulence proposed and studied 
in this work contrasts with that of iBanerjee et al.l (120071 ) . The primary reason for the 



difference is that the two works strive to answer different questions. Whereas IBanerjee et al. 
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( 120071 ) seek to determine the efficacy of isolated outflows in a quiescent media, the present 
work focuses on the efficacy of outflows toward sustaining an initial stat e of supersonic 
turbu lence. Our scenario may also be contrasted with the assumptions of iBanerjee et al. 
( 120071 ) who began their simulations with a single isothermal jet propagating into a quiescent 
ambient medium. They seek to test if a single jet can leave behind enough supersonic material 
to power turbulence on its own. Models arising from an initially turbulent state should be 
viewed as more realistic as clouds are never quiescent. Molecular clouds will be always be 
born in an environment characterized by a range of motions on different size scale due to 
formation: i.e. either gravitational collapse or collisions between large scale flows. In this 
regard we recall that both observational and theoretical studies of mole cular outflows show 
mass-velocity relations which scal e as M(v) oc v k with = 1.3 < k < 2.5 (ILada fc Fichl 11996 



Smith et all 119971 ; iLee et all 120011 ) . Thus while the bulk of material will be moving at low 
velocity much of the gas set in motion is supersonic. This is because the ambient environment 
is laced with a volume filling supersonic flow which will form a "floor" or background into 
which the outflow propagates. Thus the assumptions underlying the claims that jets can not 
drive turbulence should be considered in light of the environment s these jets i nhabit. Such 
consid eration as shown here and in the work of iNakamura fc Lil (120071 ) and ICarroll et al. 
(120081 ) shows that jets can be quite effective in driving and maintaining turbulent cascades 



We note also other factors which may cause lo wer estimates of flow speeds and mo- 
mentum imparted to entrained gas in the models of IBanerjee et all (120071 ). The outflows 



presented in this paper are driven at the high Mach number, Mj « 90 relative to the sound 
speed of the ambient gas, which is to be expected in YSO jets. The continuously driven jet 
models in this paper achieve a characteristic bow-s hock compression rati o of 150, consistent 
with the anal ytic estimate giyen by equation 5 of iBlondin et all (ll990l ). In contrast, the 
simulations of IBanerjee et al.1 (120071 ) consider isothermal hydrodynamic jets driven by flows 
with Mach numbers, Mj of 5 and 10 (with a single Mj = 20 case interacting with a cloud). 
Because the Mj = 5&10 jet models are driven with density contrast Pj/p a = 1, the inward 
and outward facing shocks at the jet Mach disk propagate at half the driving speed of the jet 
with shock Mach number M s of 2.5 and 5. Consideration of Rankine-Hugoniot jumps across 
a one-dimensional shock yields the prediction that these models would achieve compression 
ratios of M^, or 6.25 and 25, consistent with the density plots presented in that study. Thus 
lower compression ratios are achieved than is expected from higher Mach number jet models. 
Higher Mach number models exhibit more highly compressed flow structures take longer to 
decelerate to a subsonic regimes. 

Our models also consider the effects of radiative line cooling wh ich is more cha r acter- 
istic of actual protostellar outflow environments. The simulations in IBanerjee et al.1 (120071 ) 
utilize nearly isothermal flow conditions via a choice of the polytropic index 7 = 1.0001. 
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While useful for approximating the compressive effects of cooling, it is important to note the 
limitations of such an approach which can, in general, produce differences from a conditions 
in which cooling is applied. While the isothermal approximation will produce high compres- 
sions near shocks which are characteristic of flows subject to rapid radiative energy loss, it 
does not capture the effect of mechanical cooling in rarefactions. Such effects are particularly 
important for transient outflow models which employ a time-dependent jet launching. This 
effect can be illustrated by considering the one dimensional impulsive withdrawal of a piston 
with constant speed v p , from a uniform gas with sound speed c a . The leading front of the 
rarefaction has sound speed c given as, 



mm 



— z — \Vp\,0 



This simplifies to c = c Q for models which use an isothermal equation of state. While 
isothermal models adequately approximate the effect of rapid radiative energy loss and high 
compression across strong shocks, such models do not capture the effect of mechanical cooling 
across rarefactions. 

Strong rarefactions which ensue after the driving source of the outflow expire s are char- 



acteristic of outflow cavity/transient outflow models ICunningham et al.l (l2006bl ). Because 
post-rarefaction isothermal flows are not subject to the effect of mechanical cooling, the 
sound speed in these flow regions is overestimated and will likely causes a systematic un- 
derestimation of the Mach number in such regions of the flow field. This effect biases the 
results of such models in a manner that is adverse to the development supersonic motions. 
These effects should be studied in future works to ascertain their impact on the volumetric 
distribution of mach numbers in single jets driven into a quiescent medium. 



Conclusion 



The results presented in this paper show that outflow driven cavities contribute to tur- 
bulence in their immediate environment provided there exists some mechanism to disrupt the 
(mostly) laminar flow patterns that would result from an undisturbed outflow in a quiescent 
environment. The particular form of the disruption in this work is that of decaying turbu- 
lence in the ambient environment with an initial turbulent velocity that is nominally 10% 
that of the propagation speed of the outflow bow-shock. We have shown that even structures 
which are denser and propagate more quickly than the ambient turbulence are significantly 
disrupted, despite their higher momentum density. Outflo w structures can be dis rupted by 
the development of in s tabili ti es such as such as radiative ([Sutherland et al.l 120031 ) and thin 
shell IVishniac &: Ryu! (119891 ); IVishniad (119941 ) modes which have growth times far shorter 
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then the evolutionary timescales of the flow. The wide spectrum of motions present in the 
ambient turbulence provides a multi-mode spectrum of perturbation seeds which initiate the 
growth of these instabilities. Extinct outflow cavities thus become fully subsumed into the 
background turbulence in the time required for these disruptions to c ross the cavit y . Ob- 
servation of outflow driven cavities similar to the work carried out in iQuillen et al.l (120051 ) 
therefore become increasingly confused with the ambient turbulence during this transition. 
For the case of highly collimated outflow structures which are prototypical of low-mass star 
formation the relevant crossing distance is the radius of the roughly cylindrical outflow bow- 
shock. The disruption time for low mass protostellar outflows is therefore short relative to 
the turbulent crossing time of their parent cores. After the disruption time, the outflow cav- 
ities act to support turbulence in the local environment on a driving scale comparable to the 
size of the outflow cavity. We conclude therefore that protostellar outflows p rovide a efficient 
form o f dynamical feedback in low mass star forming cores as suggested by iNorman fc Silk 
jl980h . 

We postulate that the action which disrupts an outflow cavity need not be fully de- 
veloped turbulence in the ambient environment. The interaction of an extinct outflow 
cavity with another outflow or ou tflow remnant would produce similar results. The re- 
sult of ICunningham et al.l (j2006al ) showed that large radiative energy loss from the inter- 
action of multiple active outflow actually inhibits mechanical support. Here we empha- 
size cavity-cavity in teractions can also drive turbulence and this mechanism as shown by 
Carroll et al.l (120081) . Future work should also focus on ecological studies of multiple star 



formation (ILi fc Nakamural 120061 ; iNakamura fc Lil 120071 ) . Simulation of such systems over 
long timescales with sufficient resolution could determine the rate of steady-state turbulent 
decay resulting from a steady rate of outflow formation. 
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